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The usual (Bunch-Davies) Feynman propagator ol a massless field is not well defined in an expand- 
ing universe due to the presence of infrared divergences. We propose a new propagator which yields 
infrared finite answers to any correlation function. The key point is that in a de Sitter spacetime 
there is an ambiguity in the zero-mode of the propagator. This ambiguity can be used to cancel the 
apparent divergences which arise in some loop calculations in eternally (or semi-eternally) inflating 
spacetime. We refer to this process as zero-mode modification. The residual ambiguity is fixed by 
observational measurement. We also discuss the application of this method to calculations involving 
the graviton propagator. 

o : 

CN . PACS numbers: 98.80Cq 

(N ■ I. INTRODUCTION 

Improvements in cosmological observations may soon lead to the first detection of non-Gaussianity in the cosmic 
microwave background [1] or in the large scale structure of the universe [2]. These non-Gaussianities could have a 
primordial origin in theories where the inflaton has self interactions or interactions with other light fields (see [3] for 
a review) . As such there has been much recent interest in revisiting cosmological perturbation theory beyond leading 
order. In addition to non-Gaussianities, the inclusion of interactions of the inflaton naturally leads one to wonder 
about loop corrections to other observables, such as the power spectrum. 

From dimensional analysis one might expect these loop corrections to be very small, and many recent calculations 
have confirmed this expectation in several different settings [4]. Nevertheless, there is some expectation that loop 
corrections may be bigger than naive dimensional analysis would predict. The reason is that in slow roll inflation, 
the inflaton is almost massless (i.e. its mass is much lower than the Hubble scale during inflation). If the inflaton is 
, treated as a massless field, loop corrections to correlation functions are afflicted with infrared (IR) divergences due to 

■ massless particles propagating in loops, and these make the loop diagrams formally infinite as compared to the tree 
level contributions. This would lead to the breakdown of perturbation theory, and the tree level calculation would be 

■ rendered untrustworthy. A resolution of these IR divergences is essential to extract proper predictions 1 . 

These divergences have been investigated by many authors [5-9]. In principle, one may avoid the issue of these 
t-H . divergences by including deviations from exact masslessness and exact de Sitter geometries. For example, if the inflaton 
has a mass, the IR divergences will be cut off (indeed, since the mass is perturbatively renormalized, it requires a 
severe fine tuning to keep the field truly massless). It has been argued that in Aqb theory, a mass is dynamically 
generated [10-12]. The tilt in the potential in slow-roll models can also lead to a cutoff of the divergences. Finally, if 
inflation only occurs over a finite period, the boundary conditions in the far past also modify the behavior of the low 
momentum modes and cut off the divergences [13]. 

For these reasons, it might appear natural to treat the IR divergences of the massless field as unphysical, and to 
assume that one or more of these mechanisms is responsible for obtaining finite results. This is unsatisfactory for 
a few reasons. In the first place, it is perfectly possible to consider a massless scalar field that is not the inflaton, 
and which can have an arbitrary potential. There appears to be no obvious reason that the quantum field theory of 
such a scalar should break down in the exact de Sitter geometry. Secondly, even if the IR divergences are cut off by 
these physical corrections, the loop corrections may be large, and we will need to understand their effects. Finally, 
gravitons in de Sitter space also have propagators similar to those of massless scalar fields. Since their mass is not 
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1 Note that in single field inflationary models, the curvature perturbation defined on three dimensional hypersurfaces only has derivative 
interactions and loop corrections are IR finite. We do not expect this to be generally true in multi-field models of inflation or when 
computing quantities other than correlation functions of the curvature perturbations. 
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perturbatively renormalized, we are forced to deal with the divergences when we calculate corrections to observables 
which involve graviton loops. 

Here, we shall reexamine these divergences in the context of a massless scalar field theory in a (flat) Friedman- 
Robertson- Walker (FRW) space-time with constant Hubble rate (this scalar is not the inflaton) . This spacetime is a 
Poincare patch of exact de Sitter; it is locally identical to de Sitter but has only six global Killing vectors. We will 
ignore the backreaction of the scalar on the metric (implicitly assuming that Mpi — >■ oo). We argue that the resolution 
to the IR divergences involves a correct understanding of what the Feynman propagator is in an inflating spacetime. 
In fact, just as in [14], we find that a proper treatment of the zero momentum mode removes this divergence. The 
main difference between the FRW case and the de Sitter case (see [14]) is that the spectrum is continuous, and some 
care must be taken to correctly change only the zero mode. The end result is IR finite but will in general depend on 
time (which is fine for a FRW background) . We then show that this modification is compatible with other approaches 
to treating the IR divergences. However, our method is more direct and is easier for computational purposes. 

In section 2, we describe the correct treatment of the zero mode in the propagator, and in section 3 we illustrate 
with an example of a loop calculation. In section 4, we consider the graviton case. As in the scalar case, the graviton 
propagator must be modified at tree level. However, this modification applies only to the graviton zero mode, which is 
a pure gauge mode. By choosing an appropriate gauge, the IR divergences can be gauged away (as argued to be true 
in [19, 20] (see also [21, 22] for more recent discussions). This provides an IR finite method for computing graviton 
loop corrections. Our conclusions are presented in section 5. 



II. DEFINING THE FEYNMAN PROPAGATOR 



We will consider the action for a scalar in a FRW background 



S = J dtd 3 x^ (^g^dpWrf-Vtf) 



(1) 



where 

ds 2 = g^dafdx? = a 2 (r)(-dT 2 + dx 2 ) (2) 

and we take a(r) = — In the non-interacting limit with a pure cosmological constant V(cf)) = Vq, the classical 
equation of motion for the mode functions can be solved to yield 

(f) (k, T ) = -i-^- J -(l-ikr)e i ^ (3) 

where we have chosen the usual Bunch-Davies vacuum and H = ^ = J . The Feynman propagator is then 
G f {x,t;x',t') = (0\Tcf>(x,T)cj>(x',T')\0) 

When interactions arc turned on, there will be quantum corrections to this correlation function. Since 
{4>{k,T)(f){k,T 1 )) oc fc~ 3 for small k, one would expect loop corrections to include a contribution from the low- 
momentum region of the form J 4^, leading to a logarithmic IR divergence. In fact the expression (4) is infinite 
even for finite scale factor and points which are at finite separation (due to the divergence at k = 0). It has been 
suggested [10-12] that this divergence can be cured by nonperturbative effects, which can generate a mass for the 
otherwise massless scalar field. However, this non-perturbative contribution can always be canceled by a suitable 
finite counterterm, again yielding a massless particle. We can therefore assume that our renormalization conditions 
have been chosen to give us a truly massless field </>, and in this case we still must ask how this IR divergence is to be 
treated. 

Intuitively, though, it does not seem that the divergences in (4) should be physical. For practical purposes, the 
largest observable scales today exited the inflationary Hubble patch 60 e-folds before the end of inflation. The IR 
divergence in (4) arises from modes which left the horizon long before that. As these modes themselves are not 
physically distinguishable from a constant zero-mode, it would seem possible to treat virtual propagation of these 
modes in a way which removes the divergence. 

To address this problem, we note that one can add a homogeneous solution to any inhomogeneous solution to 
obtain another solution to the inhomogeneous wave-equation; these correspond to propagators with different boundary 
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conditions. The time-ordering of the Feynman propagator defines the boundary conditions: positive energy modes 
propagate forward in time, while negative energy modes propagate backward in time. However, in the massless case, 
there is a zero frequency mode. The boundary condition on this mode is not set by the causality condition, and is an 
ambiguity in the Feynman propagator. This is perhaps easier to see in position space where the propagator satisfies 
the inhomogeneous wave equation 

d^d^G f(x, t; x' , t') = -~i5 3 (x - x')S(t - r') . (5) 

This has a constant homogeneous solution, and hence the propagator can be shifted by a constant. This constant 
piece can be thought of as a shift to the 2-point correlation function, and contributes to the Feynman propagator (we 
note that since the retarded propagator is Gr{x, t; x' , r') = 9(t — t')([</>(x, t), <j)(x' , t')]), the constant shift cancels in 
the commutator, as required by causality.) 

Note that the propagator for a massless field in Minkowski space suffers from a precisely analogous ambiguity; 
there, however, the requirement from cluster decomposition that the propagator fall off at large distance eliminates 
the ambiguity. In de Sitter space (or FRW) such a requirement cannot be imposed; the propagator necessarily grows 
at large distances. 

We can therefore add a (possibly divergent) constant to the propagator to make it finite. We modify the propagator 
by putting an IR cutoff at k = /j, (and taking fi — > 0). The finite propagator is defined to be 

G%»(x, r, x>, r>) = km J ^(1 + fcVr' + ik\r - T'\)e~^ '\ e ^-*') + f (6) 

where we have added a constant / (^-^ • To cancel the leading divergence as fi — > we take / (^^f-^j = — (^) 2 m ^jr 
(here we have introduced a new comoving scale km.) The propagator in momentum space then can be written as 



G% n (k,r,r') = lim \^-(l + k 2 TT' + ik\T-T'\)e^- r '\6{k-^) + f (^K)(2Tr) 3 S 3 (k) 



(7) 



It is easy to see that the Feynman propagator is finite. Indeed, all calculations in perturbation theory will now yield 
results free of IR divergences. We refer to this procedure as zero-mode modification. Note that while the divergences 
have canceled, we have introduced a scale km which is a finite ambiguity in the propagator; this ambiguity will need 
to be fixed by an experimental measurement. A similar analysis applies for the two point correlation function. The 
usual formula for this quantity is 

f j3i, tt2 - , , 

Go(x,r,x',r') = Maf.rMaV)) = J + ikr){l - ik T ')e-^ s -^ e^^ K (8) 

The zero-mode modified correlation function is defined analogously to be 

Gg-^r.f ,r') = lim jT (g, ^(1 + i*r)(l - ikr')e-^\^-^ + f (9) 

and is also free of divergences. 

It may appear that this method of treating the zero-mode is discontinuous; one adds a constant to the propagator 
of the zero-mode, while leaving undisturbed the propagator of modes with arbitrarily small wavenumber. However, 
this discontinuity cannot be observed, as any physical experiment can only probe wavelengths smaller than the size 
of the observable universe. One cannot distinguish the manner in which small wavenumber modes are treated from 
the manner of treatment of the zero-mode. 

The effect of zero-mode modification is to effectively remove the contributions from modes below km; it therefore 
has some similarities to imposing a hard cutoff at that scale (it is important to note, however, that our scale km does 
not correspond to an actual cutoff on the mode integration. In our formalism, all modes arc integrated over, and the 
effective scale km is independent of this limiting process.) However, when we consider higher point functions, the 
integrand of the momentum integration will have additional terms, i.e. it will be of the form p-(l + aik + 02k 2 + ...). 
The IR singularity from the leading term is still canceled by zero-mode modification, but subleading terms will yield 
a finite contribution from the modes k < km. This contribution would not be present in the case of a hard-cutoff at 
km, and this feature could potentially distinguish these two approaches. 

Note that in order for eq. (7) to be a solution of eq. (5), km and fi must be comoving scales and thus independent of 
time. UV regularization on the other hand needs to be done with a UV cutoff fixed in physical scale (Ayy = fc[/ya(r)) 
and therefore the coincidence limit of the two-point function depends on time (as in [13]) 

G% u (x,t,x,t) ~ln^ (10) 
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with an arbitrary constant clq. 

This formula also offers an alternative way of understanding our modified propagator. Equation (10) implies that 
the variance of a massless field grows with time. This fact is most easily understood in stochastic inflation. In this 
approach, one treats all modes of wavelength larger than some cutoff (typically the scale of the observable universe) 
as part of a zero mode. As each mode at the scale of the horizon freezes out, it shifts the value the stochastic field 
by a fixed magnitude, but with a random phase. The value of the field thus executes a random walk, which leads to 
a variance which grows linearly with time. If inflation is taken to begin at time t — — oo, then for any finite time, we 
get an infinite variance. Alternatively [14], one can quantize the zero mode of the field separately, such that vacuum 
state is the product of a Hilbert space (for the zero mode) and a Fock space. In this formalism, the zero mode is a free 
particle. The de Sitter invariant solution is then to put the zero mode into an eigenstate of the conjugate momentum 
operator. This again leads to an infinite variance (since by the uncertainty relation, the value of the field is arbitrary.) 
This result is naturally consistent with the stochastic interpretation; in the stochastic interpretation, the random walk 
of the scalar field causes the variance to increase with time, and this can only be dc Sitter-invariant if the variance is 
always infinite. 

In the zero-mode modified propagator, this divergence is not present. Our propagator should be therefore interpreted 
as having boundary conditions where the variance is finite at finite times, but is not finite as t = — oo. Alternatively, 
we are considering a state where the wavefunction of the zero mode is not infinitely delocalized. Our choice breaks de 
Sitter invariance; by taking a limit where the constant km in (4) becomes zero, we can restore dc Sitter invariance at 
the cost of reintroducing the IR divergences. 

It is also possible to think of our zero mode modification as a limiting case of the more standard scheme of modifying 
the Bunch-Davies vacua. Indeed, as considered in [18], one way to obtain IR finite answers in a FRW spacetime is to 
modify the wave functions to 

0(fc,r) = c 1 (fc)r 1 / 2 J ff«(-fc T ) + C2 {ky/ 2 Hf^{-kT) (11) 
where H^' 2 ^ are the Hankel functions. For appropriate choices of Ci, c 2 e.g. 

ci = k-P c 2 = fk- 2 P + ^\ 1 (12) 



the divergences can be removed. The Bunch-Davies vacuum (3) corresponds to a — and c 2 = ^j. Zero-mode 
modification corresponds to the limit p — > oo, which goes over to the Bunch Davies vacuum for all nonzero k. For 
k = 0, the limit is singular; this corresponds again to the ambiguity in the zero mode. 

III. APPLICATION TO A LOOP CALCULATION 

We now apply this modified propagator to a simple loop calculation. Consider a function N(<f>) = N n + N\(j> + 
N 2 (f> 2 + .... As our notation suggest, this function could be related to the number of c- folds which in turn determines 
the scalar curvature perturbation. The two-point correlation function is 

(N(x, r)N(y, r 1 )) = N 2 + N 2 {<j>{x, r)0(y, r')> + N 2 (</> 2 (x, t)<? (y, r')) + ■ ■ ■ 

= N 2 + N 2 G (x,t;y,t')+N 2 (G (x,t;y,t')) 2 + --- (13) 

where Go is the 2-point correlator. When this is evaluated in momentum space, we find 

(N(k,r)N(k,r')) = N 2 +N 2 (^k 7 T)^k,r'))+N 2 j r)<j>{k - k' , r^k' , r')^{k - k' , r')> + ■ ■ ■ 

( d 3 k' -> 
= N 2 + N 2 G (k,T,T') + N 2 J —G (k',T,T')G (k-k',T,T') + --- (14) 

The last term can be represented as a loop diagram. If one uses the expression (8) for the correlation function, the 
loop integral will exhibit an IR divergence. In particular, if the lower limit of integration is some cutoff scale /i, then 
the divergence is manifested through the appearance of ln^, which diverges in the [i — > limit. 

We now show that by using the zero- mode modified correlation function (9), there is no IR divergence. An infrared 
divergence in the integral in (14) arises if k! <~ or k' ~ fc, in which case one factor Go is potentially divergent. The 
integral converges very quickly for k' 3> k and so the incoming momenta serve as a UV cutoff for this loop diagram. 
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Without loss of generality, consider the integrand in the limit k' ~ 0, and integrate over the modes between zero and 
k. We then find that the relevant integral is (after zero-mode modification) 

Nl -^^(k',r,r')^(k-k',r,r') ~ A 2 2 G« n (fc, r, r') / ^-G«?{k' ,t,t>) . (15) 
Jo ( 27T ) Jo ( 27r ) 

After zero-mode modification, the integral over all modes with momentum < km cancels against the zero mode, and 
we get the contribution 

A^(fc,T,r') £ ^G^r,rO~^ 2 2 Gg n (fc,r,r')(f) 2 ln^-. (16) 

We see that the contribution to the integrand from the putatively divergent regions is canceled. The integral over the 
remaining modes yields a dependence on \n(k/km), but there is no divergence. 

Note also that a shift in the quantity km, shifts the correlation function by a term proportional to Gq", and can be 
absorbed by a change in N\ . This is analogous to a renormalization group flow, where the shift in the scale leads to a 
shift in the coupling constants. In any experimental measurement, we can choose km to be the scale of observation; 
for example, in [15], one was interested in temperature correlation function in the CMB on large scales as measured 
by WMAP and Planck and hence km was chosen at a scale of order of Hq, the Hubble radius today. The ambiguity 
in the propagator is thus absorbed into coefficients like Ni, which must be measured. More generally, a smoothing 
(window) function that depends on some experimental scale must be used when comparing to any experiments; this 
smoothing procedure ends up setting the arbitrary scale km to be the observational smoothing scale 1/H in many 
calculations. Not every physical quantity needs to be smoothed; gravitational waves, for example, are not subtracted 
against any background zero mode and, as we will discuss below, our modified propagator is quite useful in this case. 

Up to now, we have been looking at a free field and modifying the propagators at tree-level. Interactions can 
introduce further infrared divergences, which will need to be canceled by modifying /(■^ La ) at higher orders. For 

example, the two point function {<j>(k, Ti)<f>(—k, r 2 )) receives a series of corrections in an interacting theory of the form 

J dndr 2 G^(k, n, r')S(fc, n, r 2 )G^(k, r', r 2 ) (17) 

where £ is a loop diagram. In general, S will contain terms of the form ln(fcjijTi) or ln(fcji{T 2 ), which will lead to 
apparent divergences when integrating over ti,t 2 . 

However, the integrand is cut off for \kr'\ 3> 1 by the oscillatory terms in Go(fc, ti, t'). This removes the divergence 
for small ti,t 2 , but leads to terms of the form ln(^S-). Such terms are finite if k > km but significantly modify the 
propagator for k <C km. To cancel the contribution from these modes, /(^ La ) will have to be adjusted order by order 
in perturbation theory. 



IV. THE GRAVITON PROPAGATOR 



The quantization of gravitons is in many respects similar to the massless scalar case. Indeed the two physical 
polarizations of the graviton satisfy the same wave equation as the massless scalar. The issues with IR divergences 
are therefore similar. 

To preserve the gauge invariance, it is preferable to work with the Lagrangian formalism and impose a gauge fixing 
condition. A suitable gauge fixing condition was found in [16]. Following them, we define 

^daPda? = i-dr 2 + dlt 2 ) 

% = K- ( 18 ) 

where 7 M „ is the unperturbed metric of de Sitter space and h^ v is the metric perturbation. Indices are raised and 
lowered with the metric "f^, at lowest order. We then impose the gauge- fixing condition 



(19) 
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and find that the equations of motion for the quadratic part of the action reduce to 

□ (^xoi) = 



□ I -x I =0 (20) 



T 



where xoo = Xoi = Xij = i>i and x = Xu + X22 + X33 + Xoo- 

The Xij therefore obey the same wave equation as massless scalars, and their propagators will be of the form (7) 
with a different constant shift to the zero- mode for each of the Xij ■ These shifts are related to the gauge symmetry; 
the gauge-fixing condition (19) imposes the constraints 

-d T Xoo + dtXm + T ~ 1 X = 
-drXoi + d 3 x^+2r- 1 x Qt = (21) 

which can be used to fix xoo and Xoi- However, there is a residual gauge symmetry which is not fixed by the above 
constraints. The residual gauge transformations can be used to generate a constant shift of the graviton perturbations 
(in fact, the IR divergence occurs because the gauge fixing term (19) does not completely fix the gauge). The correlator 
can be made finite by an appropriate gauge choice which fixes the shifts, yielding 

Gtijmn(k,T,T') = (Xij(k)Xmn(k)) 

= {9im9jn + 9in9 jm ) (jjm + ikr){\ - ikr')e ik ^-^6(k — n) — ^H 2 In ^ <5 3 (fc))(22) 

which has no IR divergence. The propagators for Xij, Xm, an d Xoo are thus free of IR divergences 2 . Furthermore, the 
ambiguity arising from the choice of km will not be present in any gauge invariant expression, since the value of km 
can be shifted by a gauge transformation. 



V. CONCLUSION 



We have constructed infrared finite propagators in inflating spacetimes by taking the usual Bunch-Davies propaga- 
tors and modifying the zero mode; a procedure we dubbed zero mode modification. We do not implicitly assume that 
all scalar fields are massive, or that inflation is not eternal. Instead, this method utilizes the fact that in de Sitter 
spacetime, there is an inherent ambiguity in boundary conditions for the zero mode of the Feynman propagator. This 
ambiguity can be chosen to cancel the divergences which usually arise from the integration of low-lying modes in a 
loop diagram. This leaves over a finite ambiguity in the propagator, which can only be set by observation. 

Physically, this modification corresponds to the choice of a state in which the variance of scalar field is finite at 
finite conformal time. In practice, the only observable effect of this choice is to change the variance (and in general, 
higher zero-mode moments) of a stochastic field, and the result is to divorce the moments of the field from the history 
of inflation before the era when observable modes froze out. The result is a coherent justification of the practical 
approach which has often been taking when dealing with loop diagrams in the quantum expansion of inflationary 
theories. At leading order, one can see that calculations performed using zero mode modification yield results similar 
to those obtained from a hard IR cutoff at some momentum km, but these approaches will differ by terms which are 
sublcading in km. Note that at a fundamental level, our message is that the description of a massless scalar field 
in dc Sitter space or its Poincare patch depends on a scale H and an extra arbitrary parameter (in the form of km 
here). 

We have moreover argued that this approach can also be used to calculate graviton loop diagrams, in which case 
the ambiguity in the propagator is in fact a gauge ambiguity. In a more detailed calculation in a theory where gravity 
is dynamical, one should be able to see that, unlike the case of the scalar propagator, observables are independent of 



2 It has been suggested [23, 24] that Gupta-Blculer quantization of the graviton also yields unambiguous finite results for observables in 
dc Sitter space. 
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this ambiguity in graviton propagator. It would be interesting to address this issue with a concrete calculation of a 
gravitational observable. 
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